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1 Introduction 



Thermodynamic Bethe ansatz (TBA) equations have been used to investigate 
thermodynamics of various kind of solvable lattice models (see for example 
lU). Traditionally, TBA equations have been derived by the string-hyposesis 
0, Several years ago, TBA equations for the supersymmetric t — J model 
and the supersymmetric extended Hubbard model were derived from the 
T-system (a system of functional relations among transfer matrices) for the 
quantum transfer matrix (QTM) p, |, ^ ||, |^, which is independent of the 
string hypothesis. Later on, TBA equations for the XXZ-mode\ in the 



regime |A| < 1 |[Tn|, the osp(l\2) model |]TT|, the osp{l\2s) model ||T2| were 



derived by the similar procedures in @]. In addition, the TBA equation 
T0[| was analytically continued |jl3[ to the one for |A| > 1. These TBA 



m 



equations contain an infinite (or finitely many) number of unknown functions, 
and thus are not always easily treated. It is significant to simplify the TBA 
equations to tractable integral equations which contain only a finite number 
of unknown functions. 

In the context of the QTM method, for models related to algebras of 
rank one or two, nonlinear integral equations (NLIE) with finite numbers 
of unknown functions were derived by Kliimper and his collaborators P, 
HI IIR ITB|, [TH] . Although their NLIE give the same free energy as TBA 



equations, the derivations need trial and error for each model, which prevents 
their extension to higher rank case. 

Another type of NLIE with only one unknown function was proposed for 
the XXZ spin chain by Takahashi |]TU| recently. One can also derive |TB 



Takahashi's NLIE from the T-system |2g, |T1 of the QTM. In view of this 
fact, we have derived NLIE with a finite number of unknown functions 



from our T-system for the osp{l\2s) model for arbitrary rank s. In 
this paper, we shall further derive NLIE for the sl{r + 1) Uimin-Sutherland 
model [24, 25] with only a finite number (the number of rank r) of unknown 
functions from the T-system [^. This is the first explicit derivation of this 
type of NLIE for a vertex model associated with sl{r + 1) for arbitrary rank 
r 

In section 2, we introduce the sl{r + 1) Uimin-Sutherland model, and 
define the QTM and the T-system |^ related to this model. In section 
3, we derive traditional TBA equations from the T-system defined in sec- 
tion 2 without using the string hypothesis. In section 4, we derive the 
NLIE ( [4.6| ), ( [4.7] ), which are our main results. The normalized fused QTM 
{Tm\v)} {a G {l,2,...,r}; m G Z>i: the fusion degree of the model) 
in the Trotter limit iV — oo defined in section 2 play the role of the un- 
known functions of these NLIE. The traditional TBA equations ( ^.11| ), ( ^.131 ), 
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( p.l6| ) are integral equations on the variables {Ym\v)}, which connects with 
the ones {Ti^\v)} for our new NLIE ( [4 .61 ) through the relation ( |3.1| ) (cf. 
Tm\v) = \imN^ooTm\v)). In this sense, we may viewed our new NLIE 
( |4.6|) as TBA equations on the variables {Tm\v)}. For the evaluation of the 
free energy, we need only the fundamental one T^^\v). Therefore adopting 
an infinite number of fused QTM {Tm\v)} as unkown functions seems to be 
superfluous. In fact we see that these NLIE ( ^.61 ) for m = 1 form a closed set 
of equations (|4.7|) , which contains only a finite number of unknown functions 
{T}"'\v)}i<a<r- For r = 1, this set of NLIE ( ^l7|) reduces to Takahashi's 
NLIE [|ig for the XXX spin chain. On the other hand, NLIE (|;6D for 
m > 2 have never been considered before, and then they are new equations 
even in the case r = 1. Using our new NLIE ( |4.7| ), we calculate the high 
temperature expansion of the free energy in section 5. Section 6 are devoted 
to concluding remarks. Many calculations in this paper are parallel with the 
ones for the osp{l\2s) case |jl2|, but we will describe calculations con- 
cisely, not so much for reader's convenience, but for basic importance of the 
Uimin- Sutherland model. 



2 T-system and QTM method 

We shall introduce the sl{r + 1) Uimin-Sutherland model |2^, and define 
the QTM g, I, 0, 1, I and the T-system m for this model. QTM analyses 
of the Uimin-Sutherland model can be found in |TB[ . The classical counter 
part of the sl{r + 1) Uimin-Sutherland model is a special case of the Perk- 
Schultz model, whose i?-matrix[2B| is given as 



Rf,{v) = v5abS^. + 5aJ^b, (2.1) 



where v G C; a,b,iy,fi G {1,2, ...,r + 1}. We define the QTM t\ 

a/3i,...,/3jv} 
'{oi,...,ajv} 



(1)/ 



((4-')!f;:::a(-)) - 



M k=i 

where = is the Trotter number and assumed to be even; 

u = —^N (-^ ^ temperature; J is a coupling constant, where J > (resp. 
J < 0) corresponds to the anti-ferromagnetic (resp. ferromagnetic) regime); 
{/Xq} are chemical potentials; the Boltzmann constant is set to 1. The free 
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energy per site is expressed in terms of the largest eigenvalue T^^^ (0) of the 
QTM (13) at = 0: 

/ = -T lim logTi^^^(O). (2.3) 

The eigenvalue formula tI^\v) of the QTM ( ^.21 ) is imbedded into the one 
Tm\v) for a fusion hierarchy of the QTM (cf. Bazhanov-Reshetikhin formula 
in H): 

a m 

T^\^)= 5^ nn^fe'^-i(«-^-2i + 2fc))' (2-4) 

where the summation is taken over dj^k G {l;2,...,r + l} such that dj^k -< 
dj-^i^k and dj^k ^ c?j,fc+i (1 ^ 2 ^ ■ ■ ■ ^ r + 1). The functions {z{a;v)} are 
defined as 



z{a; v) = ilJa{v) 



Qa.iiv-^ia + l))Qa{v-^{a~2)) 
Qa-i{v - i(a - l))Qaiv - fa) 

for a e {l,2,...,r + 1}, (2.5) 



where Qa{v) = 11^=1 ("^^^i"^); G Z>o; Qoiv) = Qr+i{v) = 1. The vacuum 
parts are given as follows 

ipa{v) = e~(p+{v + i6a,r+i)(t)-{v -i6a,i) foT a G { 1, 2, . . . , r + 1 } , (2.6) 

where (f)±{v) = (v ± iu)^ . {v]^"*} is a solution of the Bethe ansatz equation 
(BAE) 



for A; G {1,2,..., Ma} and a G {1, 2, . . . , r}. 

We can rewrite (|2.7|) as the Reshetikhin and Wiegmann's BAE in terms 
of the representation theoretical data. 

for a G {1,2,..., r}, 

where (aa|ab) = 2^a,h — 5a,b+i — 5a,b-i is the Cartan matrix of slir + 1); {6^"^} 
are the Dynkin labels which characterize the quantum space: ^^2] i ^2] 1 • • • ? 



[0,0,..., 1], = [1,0,..., 0]; {wf} are inhomogene- 

ity parameters: ^^2^-1 — ^a,iiu, — ~^a,ri{u + ^g)', g = r + 1 (the 

dual Coxeter number); (a = e t . Note that the Dynkin labels for 
the even site and the odd site are the ones for conjugate representations 
each other. The non-negative integers {Ma} should satisfy the condition: 

Ef=l bf - ELl Mbia^K) = f (<5„,1 + 6a,r) + Ma-l - 2Ma + Ma+l G Z>o, 

where Mq = Mr+i = 0. For the row-to-row transfer matrix, the above pa- 
rameters take the values: [bf \ hj'\ . . . , h^p] = [1,0,..., 0]; Wj""^ = 0; Ca = 1- 
For a G {1, 2, . . . , r} and m G Z>i, we shall normalize (|2.4| ) as Tm\v) = 



a m 



^[[[[<P-i'v ^ ^mi-v ^ V- (2.9) 

j=i k=i 

Note that the poles of Tm\v) from the functions {Qb{v)} (dress part) are 
spurious under the BAE (|2.7| ). One can show that Tm\v) satisfies the fol- 
lowing functional relation. 

T^Hv + '^)T^:Hv -'-) = xi:Uv)f:^Uv) + f^^-'\v)f^:^'\v\ (2.10) 



where 



for a G {1, 2, . . . , r} and m G Z>i, 



T^'^Hv) = 1 for a G Z>i, 



T^^+i)(t;) = e T for m G Z>i. 

Except for the vacuum part (the part made from the function (|2.6|) ), this 
equation has the same form as the sl{r + 1) T-system in [ PT[| . 



3 Traditional TBA equations 

In this section, we will transform the T-system ( P^.10| ) into the traditional 
TBA equations by the standard procedure [^, ^. Similar argument in rela- 
tion to Stokes multipliers can be found in section 5 in [BH . 
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For m G Z>i and a G {1, 2, ... , r}, we define functions: 



By using the T-system ( 2.1CI| ), we can show that satisfy the following 
F-system: 

y:.^,^^^.^. ^^SAm^-y^;^ ,,,, 

for a G {1, 2, . . . , r} and m G Z>i, 

where hb = Sa,b+i + Sa,b-i, ^o"V) = 0- 

From a numerical analysis for finite A^, u, r, we expect that a one-string 
solution (for every color) in the sector y = Mi = M2 = ■ ■ ■ = Mr of the BAE 



( |2.7| ) provides the largest eigenvalue of the QTM ( p.2[ ) at f = at least for 
the case: /ii = /i2 = ■ ■ ■ = /^r+i = (cf. Figure Hereafter we will consider 
only this one-string solution. The following conjecture will be valid for this 
one-string solution (cf. Figures ||, H). 

Conjecture 3.1. For small u (\u\ <^ 1), a ^ {1,2, ...,r} and m G Z>i, 
every zero {zm^} ofTm\v) is located near the lines Imv = ±^^^ at least for 
the case: yUi = /i2 = ■ ■ ■ = /^r+i = 0. 

Based on this conjecture, we will establish the ANZC property (Analytic, 
NonZero and Constant asymptotics in the limit 00) in some domain 

for the functions ( |3.1|) . We find that Ym\v) has a constant asymptotics: 

lim Y^\v)= ^'zfll'. , (3.3) 



where {Qm^} is given by ( p.4| ) with a formal setting Qa{v) = (l)±{v) = 1. 
{Qin} is also characterized as a solution of the Q-system p8| : 

(QL"^)' = QL'^iQLti + Qt'^Qt^'^ for a G {1, 2, . . . , r} and m G Z>i,(3.4) 

where Qo"^ = Qm'* = 1; Qm^'^^ = e r Note that Qm coincides 

with the character of m-th symmetric a-th anti-symmetric tensor represen- 
tation V{mAa) of gl{r + 1) if we set ~^ ^a, where are orthonor- 
mal basis of the dual space of the Cartan subalgebra. In particular, for 
fJ'i = = ■ ■ ■ = f^r+i = case, Qm is given by 

n(a) _ f {m + g)\m\ f {k + a){k + g - a) \'' 

^"^ \im + a)\im + g-a)\J 1 k{k + g) J ' ^'^ 
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Figure 1: Location of the roots of the BAE for s/(3) case (A'" = 12, m = —0.05, 
/^i = = = 0), which gives the largest eigenvalue of the QTM t^i\v) at 
V = 0. Both color 1 roots {v^}^} and color 2 roots {vf^} form six one-strings. 
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Figure 2: Location of zeros ofT^'iv) for s/(3) case (m = 1, 2, 3, = 12,-u = 
—0.05, /ii = = /ti3 = 0). The zeros recede from the physical strip Imv e 
[— |, |] as m increases. 
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Figure 3: Location of zeros ofTi\v) for sl{3) case (m = 1, 2, 3, iV = 12, m = 
—0.05, /ii = /Li2 = /ti3 = 0). The zeros recede from the physical strip Imv e 
[— |, |] as m increases. 
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and ( |3.3| ) has the form 



im Yi^\v) = -f 1 (3.6) 



(|3.5| ) is the dimension of V{mAa). Based on the conjecture |3]5 and ( ^73|) , 
one can see that the functions 1 + Ym\v), 1 + {Ym\v))~'^ in the domain 
Imv e [-5,5] (0 < 5 < 1) and Y}n\v) for (a,m) 7^ (1,1) in the domain 
Imf G [— I, |] (physical strip) have the ANZC property at least for the case 

/ii = ■ ■ ■ = /ir+i = 0; while Y^^\v) has zeros of order N/2 at ±«(| + tt) if 
M < (J > 0), poles of order N/2 at ±i{\ - m) if m > (J < 0) in the 
physical strip. Thus one should modify Ym\v) as 

^'i.'^H^) = ^I'^n^) jtanh ^(t; + z(l ± u)) tanh |(t; - z(l ± n))| (3.7) 
where the sign ± is equal to that of —u. By using the relation 

tanh — (f + i) tanh — (t; — i) = 1, (3.8) 
we shall modify the Ihs of the ^-system (|3.2|) as 



Y::^^.^^^^u^ '^+y^i^)^i;y^^^^^^) , (3.9) 

2 2' nfc,(i + (K!."('f))-')'- 

for m G Z>i and a G {1, 2, . . . , r}. 

This modified F-system has the ANZC property. Then one can transform 
(|3.9| ) into nonlinear integral equations by a standard procedure. 



^ogY^\v) = ±^^^^|^log|tanh^(t; + ^(l±«))tanh|(^;-z(l^ 




for a G {1, 2, . . . , r} and m G Z>i, 




(3.10) 



where Yq"'\v) = and * is the convolution. After the Trotter limit N ^ oo 
with u = we obtain a TBA equation 

^ogYi^Hv) = -^3^ + K * log |ii±^^^llil±^ ^ (^) (3.11) 
Tcoshvrt; 1 nLi(l + (^^'^ V')'" 

for a G {1, 2, . . . , r} and m G Z>i, 
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where Yq"'\v) := and the kernel is 



K(v) = \ . (3.12) 

^ ^ 2cosh7ri; ^ ^ 

One can also rewrite ( p.ll|) as 

oc 

^i^"'"*log(l + Fi"))(t;) (3.13) 

n=l 

r oo 



T cosh nv 



6=1 n=l 



where 



Jlti^) = «bSiv) - i6a,b-l + 6a,b+i)Kiv))6n,k. (3.14) 

We find that the TBA equation ( |3.13|) coincides with the one from the string 
hypothesis in ||29| after suitable modification. By using the following relations 

mm(a,d) 

1=1 

(q—a)iT 

1 sm — — — 

GJv) = - 



"^"^ - i9-a> I ^^^V. 27TV 



9 cos — — - — h cosh 

9 y 

smh — TT^k smh „ k 



(3.15) 



1=1 ^v.v.^..^2 



2 cosh I 



we can also rewrite the TBA equation ( p.ll| ) as 



+ y A-'^ * log I + ^^-'^^^ 1/-^+^^ 1 (v) (3.16) 
for a G {1, 2, . . . , r} and m G Z>i, 
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where Y^"'\v) = 0. This TBA equation (p.l6| ) does not contain {Yr 
which contrasts with (pI]),(g7[3D. We can derive the following relation from 
<^ ), dO) and ( FTSP - 



logT«(^;) = 5^G,*log(l + r/'))(. 



b=l 



r-oo sinh(A;-u) cos(A;f ) sinhf^TT^/i;) 

/ , ■ , L (3.17) 

Jo /csinh(|fc) 



After the Trotter limit N ^ oo with u = —j^, we have the free energy per 
site. 



/ = -T lim log T^^^^ (0) = J-T lim log Ti^'^(O) 



(3.18) 



where ip{z) is the digamma function 

d 



dz 



logr(z), 7 = -^(i)- 



(3.19) 



4 New nonlinear integral equations 



It was pointed out [1^ that Takahashi's NLIE for the XX Z- mo del [19 



can 



be rederived from the T-system of the QTM. In this section, we shall derive 
our new NLIE (|4]|) and (|]3) from the T-system (g30|) . 



Tm\v) has poles only at ±/3m^ = ^^^^i + iu. These poles are of 



order N/2 at most. In addition, lim|^|^oo ^m'*(f ) = Qm is a finite number 
(cf. (p.5|)). Thus we must put 



da) 
m,j 



+ 



m,j 



(4.1) 



where the coefficients 6^^, b^^j G C are given as follows. 



da) 



^ m. 



(4.2) 
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where the contour Cm'* is a counterclockwise closed loop around (3m^ which 
does not surround — the contour is a counterclockwise closed loop 
around —^m which does not surround f3m ■ Using the T-system ( |2.1(J| ), one 
can modify (|4.2|) as 



b 

m,3 



Substituting ( ^73|) into (|4.1| ), we obtain 

iV 



+ 



+ 



for a G {1,2, . . . ,r} and m G Z>i, (4.4) 

where the contour Cm"* (resp. cl^"*) is a counterclockwise closed loop around 
which does not surround — 2/3m'* (resp. 2f3m). The first term and the 
second term in the first bracket { • • ■ } in (|4.4| ) have a common singularity 
at 0; while for m = 1, this singularity from the first term disappears since 
TQ°'\y) = 1. Thus for m = 1, the contribution to the contour integral from 
the first term in the first bracket { ■ ■ ■ } in ( |4.4|) vanishes if the contour c["'^ 
does not surround the singularity at z["'^ — /^j"^ +i (cf. Conjecture ^.Ij ). This 
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is also the case with the second bracket {■ ■ ■ } in ( [4 .41 ). Therefore for m = 1, 
(|4.4| ) reduces to 



JV 

dy ^ ~ (^) ft'\y + - -,)fi''-''\y + - f ) 

dy ^ ~ (^) ft'\y - Pj"^ + ^)f[''^'\y - + t) 

/cW27rz v-y+~p'f^ T'f\y - ~P^^^ + i) 

for a e {l,2,...,r}, (4.5) 

where the contour (resp. C^""*) is a counterclockwise closed loop around 

which does not surround zf^ — jsf^ + i, — 2/3|"'* (resp. ij"-* + — i, 2(3["'^). 
Now we shall take the Trotter limit — > cxd in (|4.4| ). 



dy 1 [ Ti"),(t/ + /3i„"^-f)Tit(l/ + /3^"^-f) 
^ T^^-'\y + - ^)ri'^'\y + (3^^ - |) ] 

^/ 1 [ ri"A(y-/ji."^ + |)Tit(y-/^^"^ + i) 

^ ri°-^^(y - + f )r^°+^^(i/ - /j^"^ + f ) I 
^ r^'^^(2/-/5i"^ + ^) J 

for a G {1, 2, . . . , r} and m G Z>i, (4.6) 



where Ti^\v) := lim^^oo Tir^(^^); := lim^v^^o = r^^'-\v) = 

^ 7^^"^(f ) = 1; T^\v) = exp [ ^ ""^.^ ) ; the contour Ct^ (resp. 

C^'*) is a counterclockwise closed loop around which satisfies the condition 
\y\ < |f — f3m^\ (resp. \y\ < \v + /5m^|) and does not surround — 2/3m^ (resp. 
2Pm). In particular for m = 1, we obtain a system of NLIE, which contains 
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only a finite number of unknown functions {T^'^\v)}i<a<r 



Q 



(a) 



+ 



+ 



'c[^) 2m {v-y- P^^t}"'' {y + f3[^^ - i) 



for ae {l,2,...,r}, (4.7) 



(a) 



where 7^' 



('■+1)/ 



MlH l-Mr + 1 

e r : 



is a counter- 



the contour c'f^ (resp. C^' 
clockwise closed loop around which satisfies the condition \y\ < \v — (3^ 
(resp. \y\ < |f + l3["'^\) and does not surround z^^ — jSl""^ + i, —2i3["'^ (resp. 
z^i^ + f3["'^ — i, 2/9j"^); z["''^ = liniTv^oo ^i""*- One can calculate the free energy 



per site / by using ( |4.71 ) and the relation 

/ = j-Tiogr/^)(o). 



We need not use ([4.6|) (for m G Z>2) to obtain the free energy per site. 
In addition, {Tm\v)}i<a<r;m&>i is given in terms of the solution of ( [4.7| ) 

{T^"'\v)}i<a<r based on a Jacobi-Trudi formula pOf. However, ( [4. 61) has the- 
oretical significance since it manifests the relation between our new NLIE and 
the traditional TEA equations ([3A1| ), ( TO) , ( CT) - In fact, (U) and (pj) 
are related to the traditional TEA equations (|3.11 ), ( 3.13 ), (|3.16 ) through 
the relation (B.ll) in the Trotter limit. 



5 High temperature expansion 



In this section, we shall calculate the high temperature expansion of the free 
energy ([4.8| ) for r = 1, 2, 3 by our new NLIE ([4.7|) . As for the XXX-model 
case, the high temperature expansion of the free energy is calculated pO] by 



Takahashi's NLIE up to order 100. We assume the following expansion for 
large T: 



J. 



,n=0 



(5.1) 



In contrasts with the XXX-model case we need further assumption 

b'fHv) = > ' ^"-r..., . (5.2) 



3=0 
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where c^^j G C. Substituting into p77| ), we can obtain the coefficients 
{b^\v)}. Here we just enumerate first few of them. b^o\v) has the form 
b^Q^ (v) = log Q^i^ . Other coefficients are given as follows. 
s/(2) case: 



(1) _ 2Qi 

(1) _ 3Qf^ _ 6Qf'' 



lOQf^ 18 gf^^ 80 gf^^ 



3,2 ~ 



3Qf' 




3Q»>; 


3Q? 


22 gf^^ 


40 of' 






gf^ 


Qf 




16 gf^' 



(5.3) 



16 





2Qf 








Qf 






Jl) _ 


3Qf 




30f> 


•-2,1 ~ 


Q? 




3Qf) 


Jl) _ 
<-3,0 ~ 


10 Qf 

^in2 


18 of' 


80 Qf^^ 



-3,1 



-3,2 



22 gf^^ 40 gf^' 



+ 



(3) 
1 



gi 



(1)^ 



+ 



(1)6 



+ 



13 gf^ 

(1)3 



gf^ sgi^^" , 16 gi' 



(2)^ 



gi 



+ 



gi 
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When Q'f' are suitably chosen, (|573|) recover the known results|Q. One see 



that these coefficients are expressed in terms of the solutions of the g-system 
(|3.4| ). And thus a solution of the T-system is given in terms of solutions of 
the g-system. The T-system is a Yang-Baxterization of the g-system. Thus 
one may say that the degree of the expansion expresses the degree of a Yang- 
Baxterization. Using {hn\o)}, we can calculate the free energy (|4.8| ) and the 
specific heat C = —T^^. We have plotted the high temperature expansion 
of the specific heat for sl{3) in figure^. For large T, this agrees with the result 
from another NLIE by Fujii and Kliimper (Fig. 4 in 0]). This indicates the 
validity of our new NLIE ( [4.7| ). We note that the high temperature expansion 



for the SU{n) Heisenberg model is briefly reported in the letter ||3T[ based 
on a completely different method. 
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Figure 4: Temperature dependence of the high temperature expansion of the 
specific heat of order 8 for s/(3), J = 1, ni = h, 112 = 0, 113 = —h case {h = 0: 
smooth hne; h — 2: broken hne; h — A: dotted hne.) 
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6 Concluding remarks 



In this paper, we have derived a system of NLIE with a finite number of 
unknown functions, which describes thermodynamics of the sl{r + 1) Uimin- 
Sutherland model. This type of NLIE for s/(r + 1) of arbitrary rank r is 
derived for the first time. In particular for r = 1, our new NLIE ( ^.7] ) reduces 
to Takahashi's NLIE [|19] for the XXX spin chain. A relation between our 
new NLIE Q and the traditional TBA equations ( CT) 

which are also derived from the T-system is clarified. The high temperature 
expansion of the free energy is discussed, in which a solution of the T-system 
is given in terms of solutions of the Q-system. We expect that we can extend 



these results to other algebras by using the T-systems in [^, |3^, |3^, 

In , fugacity expansion formulae of the free energy for the XXX spin 
chain up to an infinite order is derived from the string center equation, and 
Takahashi's NLIE [T^ is rederived from this formulae. This work may be 
viewed in part as a kind of Yang-Baxterization of the results in pQf where 



power series formulae related to a formal completeness of the Bethe ansatz 
are derived from the string center equation. We can also generally recover 
40| the results in from the Q-system. Starting from our new NLIE ( [4.7|) , 
we may derive expansion formulae similar to the ones in pH] for sl{r + 1) 



of arbitrary rank r, from which an idea toward a Yang-Baxterization of the 
results in |40] may be given since our new NLIE is based on the T-system 
which is a Yang-Baxterization of the Q-system. 

There are different types of NLIE with finite numbers of unknown func- 
tions for algebras of arbitrary rank in rather different contexts El], E2[ . Their 



origin are related to Destri and de Vega's NLIE [^. To derive analogous 



NLIE from the QTM method and to study the relation to our new NLIE 
deserve investigation. 
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